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Abstract 

Generalizing earlier results on the initial data and the final fate of dust collapse, 
we study here the relevance of the initial state of a spherically symmetric matter 
cloud towards determining its end state in the course of a continuing gravitational 
collapse. It is shown that given an arbitrary regular distribution of matter at the 
initial epoch, there always exists an evolution from this initial data which would 
result either in a black hole or a naked singularity depending on the allowed choice 
of free functions available in the solution. It follows that given any initial density 
and pressure profiles for the cloud, there is a non-zero measure set of configurations 
leading either to black holes or naked singularities, subject to the usual energy 
conditions ensuring the positivity of energy density. We also characterize here wide 
new families of black hole solutions resulting from spherically symmetric collapse 
without requiring the cosmic censorship assumption. 
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1. INTRODUCTION 

The role of initial data towards determining the end state of gravitational collapse 
of a spherically symmetric dust cloud has been analyzed in some detail recently [1] . It is 
understood now, for example, that given any initial density profile for such a cloud, one 
can always choose the velocity function describing the infall of the matter shells in such a 
manner that either a black hole or a naked singularity will result as the final state for such 
a collapse. Such results, examining the relevance of the regular initial data, from which 
the collapse of the cloud commences, towards its evolution in terms of a black hole or a 
naked singularity should lead us to a better understanding of the genericity and stability 
aspects of naked singularities. This is most essential and will have important implications 
on any possible formulations of the cosmic censorship hypothesis (see e.g. [2] for a review 
of the recent developments). Even if evolutions permitting naked singularity were allowed 
by general relativity, if these were non-generic in some well-defined sense, that could lead 
to a suitable formulation for the cosmic censorship principle. 

Our purpose here is to generalize the results such as those in [1] to collapsing clouds 
with a more general form of initial data and matter. While dust could possibly be regarded 
as a good approximation to the state of matter in the final stages of collapse [3], it is an 
idealized form which does not take into account the pressures and stresses, which could play 
an important role towards determining the end state of gravitational collapse. We consider 
here general type I matter fields [4] , which include most of the physically important forms 
such as dust, perfect fluids, massless scalar fields and so on. In fact, almost all observed 
forms of matter and equations of state would fall within this general class. Having defined 
the regular initial data for the matter cloud in terms of its initial density and pressures on 
an initial spacelike surface, we examine to what possible final states the collapse of such 
a cloud could develop, while being subject to the usual energy conditions requiring the 
positivity of energy density. It is shown that given any regular distribution of matter at 
the initial epoch, there always exists an evolution from this initial data which could result 
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either into a black hole or a naked singularity, depending on the allowed choice of free 
functions available in the solution. We also consider here a model describing collapsing 
shells, which reveals some interesting features of gravitational collapse, and the structure 
of the singularity when it is visible. 

Further, our analysis here identifies wide new families of black hole solutions in terms of 
the parameters given in the initial data space or the allowed evolutions considered. Since 
the cosmic censorship hypothesis is widely useful for black hole physics, it is of crucial 
importance to isolate and identify the causes in the gravitational collapse phenomena in 
general relativity which may be responsible for the formation of naked singularities or black 
holes. The role of the initial data in this context can not be overemphasized. In fact, it 
is already known in several recent naked singularity examples that the initial data at the 
onset of the collapse does categorize the naked singular or black hole spacetimes, resulting 
from collapse. A particularly transparent example, apart from the dust collapse, is that 
of the Vaidya-Papapetrou radiation collapse models, where the rate of implosion of the 
mass (A = dm{u)/du before the formation of singularity determines whether the collapse 
would end in a naked singularity or a black hole (see e.g. [5] for a discussion). For such 
collapse models, there exist non-zero measure sets in the space of initial data, which evolve 
in either a naked singularity or a black hole. 

In other words, given the initial matter profiles in terms of the densities and pressures 
distributions, our results here indicate a specific procedure to fine tune the evolutions so 
as to necessarily produce a black hole as the end product of collapse. Such a conclusion 
no longer requires the assumption of cosmic censorship, which is turning out to be a 
diflBcult proposition to establish rigorously. From such a perspective, ours is a more specific 
characterization, in terms of the space of initial data and the possible evolutions, in order 
to generate black holes. 

2. EINSTEIN EQUATIONS AND REGULARITY CONDITIONS 

We are interested here in the problem of how a given regular initial data set evolves 
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dynamically in a spherically symmetric spacetime in the context of the occurrence and 
nature of singularities. Our main purpose is to analyze the Einstein field equations, with 
a given initial data set such as the state of matter and the velocities of the spherical shells 
at the onset of collapse for a compact object, in order to determine the possibilities of this 
configuration evolving into either a black hole or a naked singularity. We therefore consider 
the gravitational collapse of a matter cloud that evolves from a regular initial data defined 
on an initial spacelike surface. The energy-momentum tensor has a compact support on 
this initial surface where all the physical quantities such as densities and pressures are 
regular and finite. For such a cloud of sufficiently high total mass, there will not be any 
stable equilibrium configurations available, and a continual gravitational collapse must 
ensue resulting finally into a space-time singularity of infinite densities and curvatures. 

For the matter distribution of the spacetime within the cloud undergoing gravitational 
collapse, we shall not restrict ourselves to any specific form of matter but consider general 
type I matter fields. In fact, all the known and observed physical matter fields are of this 
type, except the directed radiations which are of type II. It thus follows that as far as the 
form of the matter is concerned, we have included quite general collapse scenarios which 
include all the physically reasonable matter fields, and possible equations of state. The 
type I and type II matter fields are characterized by the existence of two real orthogonal 
eignvectors, or one double-null real eignvector respectively, in a timelike invariant 2-plane 
in the spacetime. The remaining fields of type III and type IV are considered physically 
unreasonable, as they necessarily violate the energy conditions ensuring the positivity of 
the mass-energy density. Furthermore, such fields have not been observed in nature so far, 
all the observed matter fields being of type I. Thus, we would not attribute any physical 
significance or interpretation to the same presently. A type II matter distribution corre- 
sponds to zero-rest mass fields representing a directed radiation. The radiation collapse, 
corresponding to such a field can be described by a Vaidya spacetime and has been studied 
in detail to show that both the black holes and naked singularities could result in such a 
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collapse. 

Thus we concentrate here on type I matter fields, for spherically symmetric spacetimes. 
Such a matter field, in a general coordinate system, can be expressed as 

Tab \ rpa rpb , \ rpa rph . \ rpa rph . \ rpa rph j -i \ 

— AiHiiJlii + A2rj2^2 + '^3-C'3-C'3 + Mry^ry^ [1) 

where {Ei,E2,Ez,E^) is an orthonormal basis with £^4 and {Ei,E2,Ez) being timelike 
and spacelike eignvectors respectively, with A^s [i = 1,2,3,4) being the eigenvalues. For 
such a spherically symmetric matter distribution we can choose coordinates (x* = t, r, 9, 0) 
adopted to the above orthonormal frame, and the metric is written as, 

ds^ = -e^^dt' + e^^dr^ + RH^'' (2) 

where dVt^ = dO^ + sin^ Od^i^ is the line element on a two-sphere. Here z^, '0 and R are 
functions of t and r, and the stress-energy tensor given by equation (1) has only 
diagonal components in this coordinate system (i.e. we are using a comoving coordinate 
system), given by 

Tl = -p, T^^pr, Tl^pe^T^ T',^T[ = (3) 

The quantities p,Pr, and pe are the eigenvalues of and are interpreted as the density, 
radial pressure, and tangential stresses respectively. We take the matter fields to satisfy 
the weak energy condition, i.e. the energy density as measured by any local observer must 
be non-negative, and so for any timelike vector V"' we must have 

TabV^V' > (4) 

which amounts to 

P>0, p + Pr>0, P + Pe>0 (5) 

From the point of view of the dynamical evolution of the initial data at an epoch of 
time from which the collapse commences, one has a total of six arbitrary functions of r, 
namely 

iy{ti/r) = Voir), ij{U,r) ^ ^jjo{r), R{ti,r) ^ Ro{r), 
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p{ti,r) = po{r), Pr{ti,r) = pr,{r), Pe{U,r) = Pe,{r) (6) 

These functions constituting the initial data are to be specified at some initial spacelike 
surface at an initial time t = U. The dynamical evolution of such a set of initial data is 
determined by the Einstein equations, and for the metric (2) these are given by, 

^' =-" = -» ^'^="' = -1^^ 

R' 

l^'ip + Pr) ^2{pe-pr)^ -Pr (8) 

-2R' + R'^ + R^ = (9) 
G-H = l-^ (10) 

where (") and (') represent partial derivatives with respect to t and r respectively, F = 
F{t, r) is an arbitrary function of t and r, and we have put 

G = G{t, r) = e-2^(i^')^ H = H{t, r) = e'^^lf (11) 

The initial data represented by the functions i'oi'4'oi PoiProiPOo -^o ^ire not all 
independent. From equation (8) it follows that there is a relation between them as given 

by 

<{Po+Pr„) = 2{pe^ -Pro)^ -V'r^ (12) 

^^o(.)= /f^^^-^^Vr (13) 

^ ' J \ Ro{po+PrJ Po+PrJ 

Furthermore, there is a coordinate freedom left in the choice of the scaling of the coordinate 
r, which can be used to reduce the number of independent initial data to four. It follows 
that there are only four independent arbitrary functions of r constituting the initial data. 
Evolution of this data is governed by the field equations, and we have in all five equations 
with seven unknowns, namely p,Pf zv, V', -R and F, giving us freedom of choice of two 
functions. Selection of these two free functions, subject to the given initial data and the 



weak energy condition above, determines the matter distribution and the metric of the 
spacetime and thus leads to a particular evolution for the initial data. 

In spherically symmetric spacetimes, the function F{t, r) is treated as the mass func- 
tion for the cloud with F > 0. In order to preserve the regularity of the initial data at 
t = ti, we must have F{ti, 0) = 0, that is, the mass function vanishes at the center of the 
cloud. 

In this paper, we are interested only in the spacetimes which at some initial epoch t = 
ti are free from any singularities in the initial data, i.e. the evolution of the collapse must 
develop from a regular initial data. Our interest is in the evolutions of such data sets, which 
develop from this initial singularity free state, into a possible singularity of the spacetime. 
We then study the nature and structure of this singularity from the perspective of the 
cosmic censorship. The initial data is considered to be singularity free if the curvatures 
and densities are all finite, that is the Kretchmann scalar is bounded on the initial surface. 
The Kretchmann scalar K = R°'^'^'^ Rahcd for a spherically symmetric spacetime, as given 
by the metric in (2), can be put in the following form, 

K = - It^ + (14) 

where 

= (T2f = (Pr + 2po - pf (16) 

T2 = T«^T„b = p2^p2^2p2 (17) 

Thus, a singularity will appear on the initial surface if either the density p, or one of the 
pressures become unbounded at any point on the initial surface, or [F/R^) ^ oo at any 
point. We require that at the initial surface t = U, the density and pressures are finite and 
bounded. Further more, we have 

F{U,r) = J po{Ro)RldRo (18) 
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and hence the spacetime is singularity free initially in the sense that the Kretchmann 
scalar, density and pressures are finite. But as the collapse evolves, a singularity could 
develop at a later time whenever either of the density or one of the pressures become 
unbounded. We shall consider below such specific evolutions of the initial data which 
model a gravitationally collapsing matter cloud. 

3. COLLAPSING MATTER CLOUDS 

The initial data for a collapsing matter cloud basically means the initial densities and 
pressures describing the initial state of matter at the onset of collapse, namely, PoiProiPOoi 
and the function -i/^o, which is related to the initial velocity of the collapsing shells. Prom 
the point of view of the cosmic censorship hypothesis, all regular, physically reasonable 
initial data sets should evolve into a black hole. It is known, of course, that there are 
initial data sets with reasonable forms of matter such as the Lemaitre-Tolman-Bondi dust 
collapse, radiation collapse, and self-similar adiabatic perfect fiuid scenarios, which result 
into either a black hole or a naked singularity, depending on the nature of initial data and 
initial distributions of densities and pressures. Thus, the above broad concept of cosmic 
censorship has to be fine tuned and made more precise in order to arrive at a genuine 
formulation and proof for censorship (see e.g. [2] or [6]). 

It follows that either for such a purpose of cosmic censorship, or to understand the 
naked singularities of gravitational collapse better, we need to examine the evolutions of 
an arbitrary, but regular and physically reasonable set of initial data undergoing the grav- 
itational collapse. This requires defining more precisely as to what constitutes a regular, 
physically reasonable, initial data. As far as the form of the matter is concerned, since all 
the observed forms of matter have been of type I as discussed earlier, we have considered 
here general type I matter fields. Next, it is widely accepted that all physically reasonable 
matter forms must satisfy the energy conditions ensuring the positivity of mass-energy 
densities. Therefore, one of the conditions we must impose naturally is that at the onset 
of the collapse, all the initial data sets specifying the density and pressures profiles of the 
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cloud must satisfy an energy condition, and the same must hold during the later evolution 
of collapse. The regularity of the initial data also means that at the onset of collapse at 
some t — ti, the spacetime must be singularity free, as discussed earlier. Further more, 
from the point of view of general relativity the data should be at least differentiable. 
Therefore we require that the functions Po^Pr^^Peo be atleast function of r. The same 
applies to the metric functions and '^o which must be functions of r throughout the 
spacelike section t = ti for all r > 0. 

As stated above, for the initial data to be physically reasonable the weak energy 
condition is to be satisfied. This imposes the following restrictions on the initial data, 

Po{r) > 0, Po+Pr„ > 0, Po{r)+pe„{r) > (19) 

Further, for the evolution of the matter in the spacetime to be physically realistic, energy 
conditions (5) have to be satisfied. Therefore we assume that through out the cloud, at 
the initial and at all later epochs the weak energy condition is satisfied during the entire 
evolution of the collapsing cloud till the formation of the spacetime singularity at R = 0. 

Since the function Vo is related to the initial matter data by equation (13), and is 
a function of r, this puts certain restrictions on the arbitrariness of the choice of the 
functions PoiPvoiPOo- -^^^ example, we have 

be„-PrJr=o = (20) 

For the sake of physical reasonableness, we require the center r = to be the regular center 
for the cloud, which means Rit, 0) = 0. Also, one would like to have the initial density 
Po(0) > at the center r = 0. This implies that we have 

Po(0)+p.„(0)>0, Voir) = r''h{r) (21) 

where h{r) is at least a function of r for r = 0, and at least a function for r > 0. 
This means that the pressure gradients vanish at the center r = 0, basically meaning that 
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the forces vanish at the center. In this section we consider the above scenario for the 
sake of physical reasonableness, however, it is possible to give a more general formalism 
independent of requirements such as above, as we shall discuss in section 6. 

In fact, it would be reasonable to require that pressures be positive at the onset of 
the collapse, since for astrophysical bodies physically we would prefer the pressures rather 
than tensions. Further more, to make the scenario physically more appealing, we may 
require the density to be decreasing as we move away from the center r = 0. In that case, 
for various reasonable equations of state such as p = kp,0 < < 1, (a perfect fluid), or 
p = kp'^ , the pressure also will decrease away from the center together with the decreasing 
density. This will be typically the case in the massive bodies such as stars and such other 
astrophysical systems. Then as such the energy conditions could impose restrictions on the 
maximum size of the matter cloud with such an initial density and pressure distribution. 

We use the coordinate freedom available in rescaling the radial coordinate r such that 



The physical area radius R then monotonically increases with the coordinate r, and there 
are no shell-crossings on the initial surface, with R' — 1. Since we are considering gravita- 



Consider now the gravitational collapse of a matter cloud with a general initial data 
as prescribed above. We supply the two free functions F and u as below. 



From equations (7) to (13), the evolution of the collapse is described by the equations. 



R{U,r) = Ro{r)=r 



(22) 



tional collapse, we also have R < 0. 



V — 



c{t) + iyo{R), F^f{r) + F,{R) 



(23) 




(24) 



G = b{r)e 



(25) 



Vrr 



= —a 



(t)e^° ^Jh{r)Re''o -R + f + Fo 



(26) 
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Pr = (28) 
2p0 = Ru,R (p + Pr) + 2pr + Rpr,R (29) 

Fo{R) = -l r^PrJr = -R^p{R), f{r)^f {po + Prjdr = rh{r) + r^r) (30) 
./o Jo 

Here F{ti, r) — 2r^e{r). The quantities e(r) and p{r) are to be treated as the average 
mass and pressure densities of the cloud, and are decreasing functions of r. Since pr„ is a 
positive function on the initial surface, it follows that < throughout the spacetime, 
and as such the radial pressure is non-negative throughout the spacetime. The arbitrary 
function h{r) characterizes the velocity of the spherical shells at the initial time t = U. 
We are dealing with the collapse situation with ^ < 0, therefore the arbitrary function 
a{t) > 0. 

It should be noted that for the case Pr = Pe = 0, the above set of equations, and 
the collapse models reduces to the Tolman-Bondi-Lemaitre case of general inhomogeneous 
dust collapse. In that case, the results in [1] apply as far as the role of the initial data 
towards determining the final fate of collapse is concerned in terms of either a black hole 
or a naked singularity. Thus, the models here may also be viewed as directly generalizing 
the Tolman-Bondi-Lemaitre results to include both the radial and tangential pressures in 
order to investigate again the role of initial data towards the final fate of collapse. Certain 
examples have been examined which include non-zero pressures, and it is seen that again 
both black holes and naked singularities arise as final fate of collapse [7]. 

Since r = is the regular center of the cloud meaning R{t, 0) = 0, it follows from 
equation (26) that 



V^v = -a{t)e''°^/v^{h{R)b{r) - p{R)) + bo{r)v + e(r) + p{r) (31) 

where the arbitrary function b{r) — 1 + r'^bo{r), such that bo{r) is at least a function 
of r for r = 0, and a function for r > 0, and we have put 

R^rv{t,r), v{U,r) = l (32) 
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h{R) = h{rv) = -5-5 (33) 

The functions bo{r),h{rv),v{t,r),u{rv), and fo{r) are all at least functions of their 
arguments. Note that ai t = U we have v = 1 and since v < 0, we have v < 1 throughout 
the spacetime. 

We note that the quantity R{t, r) > here is the area radius in the sense that 
4:TTR^{t,r) gives the proper area of the mass shells at any given value of the comoving 
coordinate r for a given epoch of time. The area of such a shell at r =const. goes to zero 
when R{t,r) = 0. In this sense, the curve t = ts{r) such that R{ts,r) — describes the 
singularity in the spacetime where the mass shells are collapsing to a vanishing volume, 
where the density and pressures diverge. This shell-focusing singularity occurs along the 
curve t — ts{r) such that v{ts,r) = 0, the Kretchmann scalar diverges at such points. 
Using the remaining degree of freedom left in the scaling of the time coordinate t we could 
put a{t) = 1. Equation (31) can then be integrated with the initial condition v{ti,r) = 1 
to obtain the function v(t,r). Note that the coordinate r is treated as a constant in the 
equation. We in fact get 

- 1 (34) 



V \/ho{r)ve^'^° + e^^° {v^{h{rv) — p{rv)) + e(r) + p{r)) 

where we have chosen for the sake of simplicity U = 0. The time t = ts{r) corresponds to 
the occurrence of singularity is then given by, 

/ X \fvdv 
t = ts{r)= . . ... . ... . , . . , = (35) 



\/ho{r)ve^^° + e^^° {v^{h{rv) — p{rv)) + e(r) + p{r)) 

For such a collapsing cloud, a singularity can also occur at points where i?' = 0, which 
are termed the shell-crossing singularities. Though these are singularities of a weaker na- 
ture in general through which the spacetime can possibly be extended using a suitable ex- 
tension procedure [8] , the comoving coordinate system we have used here may break down 
and the metric might possibly become degenerate at the points where R' — 0. However, 
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our purpose here is to study the shell-focusing singularity at i? = 0, which is essentially dif- 
ferent and could be much stronger gravitationally as compared to the shell-crossings which 
are possibly delta-function like singularities, caused by different shells crossing each other 
where the density momentarily blows up. Hence, we choose the evolution of the initial 
data in such a manner that any shell-crossings are avoided in the collapse, except possibly 
at the singularity. Here it may be useful to mention that a similar situation regarding the 
occurrence of shell-crossings arises in Tolman-Bondi-Lemaitre dust collapse models also. 
However, as has been pointed out in earlier works, for a given initial density profile one can 
always choose appropriate initial velocity of the dust shells such that during the evolution 
no shell-crossings are encountered or visa-versa. 

In a similar fashion, this can be achieved in the present general case as well by a 
suitable choice of the functions involved, as specified below. At a given epoch of time, the 
functions e(r), p(r), Uo{R) and h{R) are at least functions and further more e(r) and 
p{r) are decreasing functions of r. The function bo{r) is an arbitrary function representing 
the initial velocities of the collapsing spherical shells. From equation (35) it is clear that 
the singularity time ts{r) is an explicit function of the velocity function bo{r) which is a free 
function, and one can choose it in such a way that ts{r) is an increasing function of r, i.e. 
dtg/dr > 0. The exact nature of such velocity functions bo{r) for which ts{r) in equation 
(35) is an increasing function of r, depends upon the exact behavior of initial density and 
pressures within the cloud. For example, for a matter cloud initially satisfying an equation 
of state of the type p = ap^ one of the many possibilities is the function ho{r) > such 
that b'g{r) < and is less than a certain minimum for r(, > r > 0. For all such functions 
bo{r), therefore, the singularity curve t = ts{r) is an increasing curve for all allowed values 
of coordinate r, and the successive spherical shells within the cloud collapse to singularity 
successively, and shell crossings do not occur. Thus R' = v + rv' > (note that R' = 1 
initially) throughout the spacetime. Furthermore, if [r^6o]' > than ^/vR' < 1 within the 
cloud for 1 > f > 0. 
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The shell-focusing singularity R = occurs first at r = and the time of occurrence 
of such a singularity from equation (35) is given by, 

ts. = ts{0) = (36) 
Jo y/v^{ho-po) + booV + eo+Po 

where ho = ho{0)-iPo — p(0),eo = e(0),6oo = ^o(O) are constants related to the central 
density and pressures, In fact, near the center r = we have 

ts{r)^t,^+rX{0) + O{r^)... (37) 

where the function X = X{v) is given by 



x_ i ^o+Po + bpoV + {hp - po) Vvjei +pi+ biv - v'^hi)dv 

eo+Po A {vHK-po)+booV + eo+Po)y^ ^ ^ 

where ei = — e'(0),pi = —p'{0),bi = — 6^(0),/ii = h,ji{0). For the case that we are 
considering where the pressures have been taken to be positive, the central singularity 
at r could be naked, all subsequent singularities with r > are covered as the 
quantity F/R ^ oo and the trapped surfaces and the apparent horizon develop prior to 
the formation of the singularity. (When the pressures are allowed to be negative, still 
subject to the validity of the weak energy condition, the other parts of the singularity can 
be visible in principle, see e.g. Cooperstock et al in [7]). It thus remains only to examine 
the nature of the central singularity. 

We now show that the weak energy condition is satisfied by the class of collapse models 
under consideration. For the models to be physically reasonable, it is necessary that the 
positivity of energy density is preserved at the initial epoch, and also throughout the entire 
evolution of the collapse, till the formation of the spacetime singularity at i? = 0. At the 
initial epoch, for the regularity and physical reasonableness of the initial data, we require 
the initial matter density Poir) to be positive, and also the initial radial and tangential 
pressures Pr^ and pg^ to be positive, or at least non-negative. Thus we have, 

Po{r) > 0, pr^ > 0, pe^ > 
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The initial density Po(0) = Pc is maximum at the regular center r = 0, and decreasing 
away from the center, and the central pressures Pr,Pe are also positive. It is therefore 
intuitively clarifying to understand the behavior of these physically relevant quantities in 
the neighborhood of the regular center r = 0, and near the singularity, during the evolution 
of the collapsing cloud. Firstly, it is clear from equation (28) that the radial pressure pr 
will be positive throughout out the subsequent evolution of the cloud. It then follows from 
equation (27) and (28) that p + Pr > throughout the spacetime, since it is that way to 
begin with. Next, note that in the near regions of the center at r = 0, and that of the 
singularity at f = 0, we have p > 0, p + pr > 0,p + pe > 0, if the same is satisfied at the 
initial surface, regardless of any restrictions on the form of initial densities and pressures. 
Actually, the initial density Po(0) = pc is maximum at the regular center r = 0, and 
decreasing away from the center, and the central pressures Pr,Pe sixe also positive there. 
In the near regions of the center r = 0, ^/vR' = v^/^ + rX{v) + O(r^), X{v) < 1, and we 
have from equation (27) to (29) 

v^R'p = pc- Pir + Pr^l -v^- v^/^X{v)r) + 0{r'^) (39) 

v^R'ip + pe) ^pc+ pe. - Pir + 0{r^) (40) 

Pc + Pr. 

where pci Pr. and pg^ represent the central density and pressures at the initial epoch 
t = ti, and pi is the first derivative of Po{r) at the center. As v we have p +oo, 
p + Pr ^ +00, and p + pe ^ +oo. Thus p > 0, p + Pr > and p + pe > 0, and the 
energy conditions are clearly satisfied in a certain neighborhood of the regular center r = 
through out the evolution as seen from the above. It is seen that at the shell focusing 
singularity v = the density becomes infinite where the curvature singularity occurs. 
In the case when the densities and pressures are initially decreasing functions of r, then 
throughout the cloud they have maximum values at the center, and are decreasing within 
the cloud. 

In fact, a sufficient condition that the energy conditions be satisfied for a matter cloud 
of a finite size Tc > r > at all times is that (/(r) — f{R))' > within the coordinate range 
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Tc > r > 0. Note that the inequahty is satisfied both initially for v = 1 (i.e. at t = 0), 
in near central region throughout the evolution, and in the neighborhood of the entire 
singularity curve v = 0. In fact, as v ^ we have p — > +oo, p + p^. — > +oo, p +P0 ^ +oo. 
Further, it follows from equation (30) that for any given value of the radial coordinate r, 

pr pR 

f{r)= r\po+Pr)dr>f{R{r,t))= r\po+Pr)dr 
Jo Jo 

Hence the mass function F > for all r at all times throughout the cloud. Since /' = 
r^(po + Pro) > and is an increasing function of r for Tc > r > where f"{rc) = 0, it 
follows that f'{r) > f^ii{R)R', and energy conditions are satisfied within rc > r > 0. Thus, 
for a given initial matter distribution as above, there is a finite cloud where the energy 
conditions are satisfied through out the evolution of the initial data. Next, for a cloud 
of arbitrary size and for situations with completely arbitrary initial density and pressures 
which are not necessarily decreasing (or even including negative pressures, or pressures 
and density increasing with increasing r etc.), note that the negative contribution in the 
expressions for p and p + pe is finite throughout the evolution while the positive terms are 
maximum at the center and become unbounded at the singularity. The behavior of the 
term R'f,R — r^v^R'{p(^R) + pr{R)) (specificely R' oc b~^{r) for v > 0) depends explicitly 
on the velocity function bo{r) which is a free function, while f'{r) = r'^{po{r) + Pr{r)) is 
independent of the same. Therefore a suitable choice of bo{r) such that 

Prir)+Pr{r) 



Po{R)+Pr{R) 



>v^{v + rv') (41) 



within the range rb > r > would result in the energy conditions being satisfied for such 
an arbitrary size cloud for a given initial data. 

Clearly, depending upon the conditions and values on the boundary of the cloud for 
the density and pressures, the exterior matching for the general class of models given here 
could be Schwarzschild, or one of the other suitable spacetimes, such as for example a 
Vaidya model spacetime. We shall not enter into these details of the matching presently, 
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as our main concern here is the local nakedness or otherwise of the central singularity at 
i? = 0. 



4. APPARENT HORIZON AND THE NATURE OF SINGULARITY 

The apparent horizon within the collapsing cloud is given by R/F = 1, which gives 
the boundary of the trapped surface region in the spacetime. It is the behavior of the 
apparent horizon curve (which meets the central singularity ai R — r — 0) near the center 
which essentially determines the visibility, or otherwise, of the central singularity. For 
example, it is known within the context of the Tolman-Bondi-Lemaitre models that the 
apparent horizon can be either past pointing timelike or null, or it can be spacelike, as can 
be seen by examining the nature of the induced metric on this surface. This is unlike the 
event horizon curve which is always future pointing null. If the neighborhood of the center 
gets trapped earlier than the singularity, then it is covered, and if that is not the case the 
singularity can be naked [9], with families of nonspacelike trajectories escaping from it. 

In order to consider the possibility of existence of such families, and to examine the 
nature of the central singularity occurring at i? = 0, r = in the general class of models 
considered in the previous section, let us consider the equation of the outgoing radial null 
geodesies which is given by, 

^ = e^-" (42) 
dr 

The singularity appears at the point vitg^r) = which corresponds to R{ts,r) — 0, 
therefore if there are future directed outgoing radial null geodesies, terminating in the past 
at the singularity, then along these trajectories we have — > as t — > t^- 

Writing the equation for these radial null geodesies in terms of the variables {u — 
r^/^, R) we get, 

dR 3 , R ^,1, 1 — /.„s 
— = -- + v^^V ^ ^ (43) 

If the null geodesies terminate in the past at the singularity with a definite tangent, then at 
the singularity the tangent to the geodesies dR/du > in the {u, R) plane, and must have 
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a finite value. In the case of collapsing ball of matter we are considering, all singularities at 
r > are covered since F/R —>■ oo, and therefore dR/du —>■ — oo, and only the singularity 
at center r = could be naked. As mentioned earlier, for the case when R' > near the 
central singularity, we have 

Xo= lim -= lim ^^xl/' = -X{0) (44) 

where -^(0) is given by equation (37). Since -^(0) > the singularity is at least locally 
naked. The behavior of outgoing radial null geodesies in the neighborhood of the singularity 
are described hy R — XqU in {R, u) plane and in (t, r) plane it is given by 

t-t,(0)=xor^/^ (45) 

The global visibility of such a singularity will depend on the overall behavior of the 
various functions concerned within the matter cloud and we shall not go into those details 
presently. It has been seen, however, from the study of various examples so far, that once 
the singularity is locally naked, one can always make it globally visible by a suitable allowed 
choice of functions. Note that in cases where the choice of ho{r) is such that X(0) < the 
singularity would be covered. When X(0) = 0, one has to consider the next higher order 
expansion term which is nonvanishing in equation (37), and that will then determine the 
nature of the singularity by means of essentially a similar analysis. 

One can also write the equation for these radial null geodesies in terms of the variables 
(t, R) to see how the area radius R grows along these outgoing null geodesies with increasing 
values of time. As mentioned earlier, for the case where i?' > near the central singularity, 
we get 

Xo = lim ^Vt = lim ^ == hm [e" ^ = 1 (46) 

t^to,r-*0 t - ts{0) t^to,r^O dt t-^to,r^O ^ + ^/H 

Again this shows that the singularity is naked at least locally. In fact, as pointed out 
above, it follows that for 6^(0) ^ the area coordinate behaves as i? = const. r^^"^ near the 
singularity. 
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5. COLLAPSING SHELLS 

We next consider here collapsing thick shells of matter, subject to the weak energy 
condition, but otherwise with a general form of matter, as this reveals some interesting 
features regarding the structure of the singularity. We rescale the coordinate r, contrary 
to the earlier case of collapsing matter clouds, such that at some initial t — ti the density 
distribution of shells is expressed as 

Po(r) = pe(l-piO (47) 

where pc > 0, pi > and a> 1 are constants. It is seen that in the scaling given as above, 
the area radius R again increases with the increasing coordinate r. 
We choose functions F and G such that 

F = f{t)R, G = a(r)e*(^) (48) 

where f{t) and a(r) are arbitrary functions of t and r only. Prom field equations (7) 
to (11) we get the matter distribution in spacetime being characterized by p,Pr and pe 

e® = ae* + 1 - / (52) 
R = -b{t)e^ Vae* + 1 - / (53) 

where H = and Q{R,t) is treated as function of R and t and represents partial 
derivative with respect to R with t kept constant. Also, b{t) > is an arbitrary function 
of t. If the spacetime characterized by the solutions of the field equations given by these 
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equations does indeed evolve from the initial matter data po{r),Pr^{r) aindpg^{r), one must 
have at the initial epoch t = ti 

p{t,,r) ^ Po{r) = fo^fiU), (54) 

Pr{ti,r) =pro = fi = [f]t^ti,Vo{r) = -v{ti,r), (55) 

Pe{ti,r) = 77^(1 - RofJ^,R„ ), MRo) = Q{Ro,ti)-logCo (56) 
where Cq is a constant. On further simplification we get 

Pc - Pir"' 



^^^|^e.e-M-^), mo) = ^o (58) 

a(r) = (c,e2^(^°) - 1 + /,)e-*° (59) 

where e{Ro) = Po+Pro and b{ti) = bo, Cq, fo, fi are constants. Free constant Cq > is such 
that a(r) > for all allowed r. The functional form of ^{R) is given by 

Note that ^{R) is a function of R for all i? > and o(r) is also a function of r 
within the cloud. These conditions further imply that e® is also a function of R and t 
for > 0. Differentiating (53) with respect to r we get 

(ee),H-(e^),K ^ ^ 

At t = i?^ > 0, and so from the above equation it follows that R' could vanish only when 
either e* = or at points where e — » oo. Since both e"^ > and > are functions 
for i? > 0, it follows that R' could vanish only at i? = 0. Thus for t < tg, we have R' > 0. 

The inner shell of the matter is labeled by r = and the outermost shell by r = r^. At 
t = ti the density and pressures are finite throughout the cloud and the Kretchmann scalar 
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is bounded at t = tj. The singularity appears first when the innermost shell collapses to 
zero volume in the sense that -R(ts, 0) = 0. The time of collapse of further successive shells 
to singularity is given by R{ts{r),r) = 0. The weak energy condition is satisfied if f{t)>0 
and f{t) > 0. The interesting feature of such an example is that at the singularity R = 
the mass function F = 0, and as such for f{t)<l we have R> F and the entire singularity 
curve is not covered by the apparent horizon. Thus we have a naked singularity where not 
just the central point but the entire singularity curve is visible. On the other hand, for 
f{t) > 1 the singularity is covered and we have a black hole. Of course, in either case, 
the mass function F vanishes on the entire singularity curve and so this is a "massless" 
singularity in that sense, but what is important is this example illustrates the possibilities 
inherent within the Einstein field equations towards the determination of the final fate of 
gravitational collapse of a massive matter cloud. 

6. CAUCHY PROBLEM AND THE INITIAL DATA 

In previous sections different classes of evolutions for the initial data leading to ei- 
ther a black hole or a naked singularity were discussed. We shall make now some general 
comments on all possible dynamical evolutions of a given regular initial data in the con- 
text of the occurrence of the singularities, naked or otherwise, for general matter fields. 
This is related in a way to the Cauchy problem in the context of general relativity. The 
difference, however, is that while most discussions on Cauchy problem in relativity avoid 
singularities, our purpose is precisely to examine the possible evolutions of an initial data 
into a spacetime singularity, in terms of its visibility or otherwise. 

We pointed out in [10] that a naked singularity would occur generically in spherically 
symmetric collapse situations. The point of emphasis here would be to examine this issue 
in the context of a given arbitrary initial data. We examine if the in built freedom of choice 
in selection of an equation of state in general relativity is sufficient enough to allow generic 
evolutions of an arbitrary initial data into a naked singularity or a black hole subject to 
the energy conditions being satisfied throughout the evolution. The basic issue is, firstly, 
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given an arbitrary physically reasonable initial data describing the state of the matter at 
the onset of gravitational collapse, does general relativity generically allow the evolution 
of this data into a singularity, subject to the energy conditions being satisfied throughout, 
and if so, do all these evolutions have to necessarily end either as black hole or a naked 
singularity, or is there still some degree of freedom available in the field equations which 
would allow the formation of either as desired. 

To address this issue along the lines of [10], the field equations given in equations (7) 
to (11) can be put in the following form if one uses R and t as variables instead of t and r. 



koR^ ^' koR^ R'^T 



Q,R \ 



2 



^^~F„Q{F^RR+2peR) ^^^^ 



QHl-i) + F,-, 
F,t{F,RR+2p0R) 



Q,RR-Q,t = fi{R,t,Q,Q,R) (65) 



H = G-l + ^ (66) 

where ,r and ,t represent partial derivatives with respect to R and t keeping t and R as 
constants respectively, and we have put 

Q= , =. p = F,KR+2peR (67) 

G-l+l 

MR.t, g, Q.r)-q^ [Q:r +Q + ^^ + ^j + 2^"^ ^^^^ 

The mass function F and the tangential pressure po are free functions of R and t, 
and give the distribution of matter throughout the spacctinic. As usual R at the initial 
t = is scaled as i?(0,r) = r using the scaling freedom available. Knowing F{R,t) 
and po{R,t) one integrates the second order parabolic partial differential equation (65) 
with appropriate initial and boundary data to obtain Q{R,t). Then e^^ is immediate 
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from equation (63). Equation (64) is then integrated for R{t,r) with the initial condition 
-R(0, r) = r. Equations (66) is then solved for the remaining unknown ip. The initial 
data, as pointed out in section 2, consists of four independent functions of coordinate r 
at an initial spacelike slice t = when the collapse commences. These are basically the 
density distribution Po{r), radial pressure Proi''^): and tangential pressure peo{r), which 
describe the state of matter at the onset of collapse, and the velocity distribution function 
Vo = Vo{r) = — -R(0, r), which describes the radial velocities of the spherical shells towards 
the center. Thus, for an arbitrary initial data at t = we have 



F{R, 0) = F(r, 0) = Fo{r) - / r^p^dr, F,^ {R, 0) = Fi(i?) = Fi(r) = r^ProVo (69) 



The functions F and pg are free functions to be chosen, subject to their initial values 
at t — 0. The functional behavior of F and pe for t > is completely free except that 
they should at least be C^, and such that energy conditions are satisfied throughout the 
collapse, i.e. p > 0-, p + Pi > 0. This will be the case when 



Thus energy conditions are satisfied for all the evolutions where the mass function is a 
monotone increasing function of both R and t, with non-negative tangential pressures. 

We therefore consider only those sets of arbitrary functions F{R,t) and p0{R,t) for 
which F > 0, FjR > 0, F,t > 0, pg > — p, and p = F^fm +2pgR > everywhere except per- 
haps at the singularity R = where it could vanish. Note that in case p = the parabolic 
partial differential equation (65) is replaced by an second order ordinary differential equa- 
tion for T — R with Q = c{t)F,t, and all our conclusions here would again apply in this 
case. We would further consider only those sets of functions for which /i in equation (68) 
is at least C° function of its arguments. From equation (14) which gives the expression for 
the Kretchmann scalar, it is clear that at the initial time t = it is finite since the initial 




P0{R, 0) = p0{r, 0) = peo{r) 



(70) 



F,R>0, F,t>0, pg> 



F\r_ 
koR^ 



(71) 
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data in the form of density, pressures, and vefocity functions are all at least functions, 
and therefore at t = we have F{R, 0) = R^h{R), where h{R) is at least function of 
R. Therefore, for all sets F{R,t) a shell-focusing singularity R{ts,r) = would develop 
at a later time t = ts{r) if F ~ i?", n < 3. Similar situations would occur when we have 
R'{t, r) = 0, and a shell-crossing singularity would occur. Thus for a given set of arbitrary 
initial data of matter and velocity functions at the initial spacelike hypersurface, there are 
infinite many evolutions satisfying the energy conditions, characterized by the free func- 
tions F and pg, resulting into a singularity. The question that needs to be answered is 
for a particular evolution determined by a particular choice of functions F and pe , would 
there be a naked singularity or a black hole. 

Note that the second order partial differential equation above for Q is a quasi-linear 
parabolic equation in variables R and t. The coeflBcient of Q,rr is positive and at least 

(depending upon the choice of F and pg), while that of Q,t is —1. Furthermore, the 
driving term fi{Q, R,t,Q,ji) is also at least C° function of its argument. The initial 
value of Q{R,0) = Qo{R) at t = is fixed by the initial data as mentioned earlier. 
Thus we have a well posed Cauchy problem with the boundary condition at some t — 
still free to choose. The quasi-linear parabolic partial difi'erential equations, such as (74), 
have been studied quite extensively for existence of solutions for an initial condition at 
t — 0,Q{R,0) — Qo{R), and the boundary condition at some t — ts, and the solutions 
do exist with quite general form of boundary conditions Q{R,t) = g{R,t). Therefore, for 
all sets of functions {F,pg) for which the coefficient of Q,rr in equation (65) is positive, 
and together with /i is at least C° function of its argument, the solution with quite 
general boundary conditions would exist which would not only evolve from a given initial 
data but also satisfy the energy conditions throughout evolution. This is the case as 
shown in examples considered in earlier sections. The point is that with the choice of a 
particular initial data set, namely {po,Pro,Peo, Vo), and a particular set of evolutions given 
by {F{R,t),pg{R,t) we still have the choice of the boundary condition aX t = ts at our 
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disposal. The freedom left in the form of the choice of Q(-R, tg), once chosen, completes the 
solution of the field equation and it is this choice which effects the nature of the singularity 
occurring at t — tg, R = in terms of being naked or otherwise. This is precisely what 
happens in the classes considered in earlier sections where an evolution from an arbitrary 
initial matter data can form either a black hole or a naked singularity. 

To clarify this, let us briefly consider the radial null geodesies equation, which is, 

dr dt ^ ^- ~ ^ ' ' ^ ' 



where X = R/t — tg- Since at the initial t = 0, R = r and R' = 1, and ^ < throughout 
the spacetime, the singularity would occur at a later time t = tg > 0. Prom equation 
(74) it follows that for a choice of functions {F^po) such that F > 0,p > throughout 
the spacetime, Q has at least two continuous i?-derivatives and that p has at least one 
continuous R-derivative, then i?' > and is finite throughout the spacetime if it starts 
that way, except possibly at the singularity i? = 0. Thus in such cases shell-crossings 
would not occur, and the first shell-focusing singularity would occur at the center at 
r = 0,t = tg. Note that in case F(Q,ts) > the singularity would be covered. The 
singularity would be naked and there would be outgoing radial null geodesies if the equation 
V{X) = U{X, ts)-X ^ 0, where 

Xo = lim ( = lim ( ^) = U(Xo, tg) (73) 

(R^O,t^ts) \t — tg J {R^O,t^ts) \ dt J 

has a real positive root X — Xo- In the case when there are no real positive roots for 
V{X) = 0, then the singularity is covered. The value of the function Q{R,tg) which 
appears in the root equation above thus becomes important. For all functions Q{R,tg) 
chosen such that the equation has a root gives rise to the naked singularity, the choice 
otherwise leads to a black hole. Q{R, tg) is precisely the boundary condition in equation 
(65) which one has a freedom to choose for obtaining a particular solution. Therefore, 
as illustrated by the two examples in the preceding sections, there would be a non-zero 
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measure set of evolutions from an arbitrary initial data specified at the onset of the collapse, 
which could always develop into a singularity, naked or covered. 

7. CONCLUDING REMARKS 

We considered here various classes of collapsing evolutions from arbitrary initial data 
and it is shown that these can be chosen so as to form either a black hole or a naked 
singularity. Such an initial data, prescribed in terms of the initial density and pressure 
profiles of the cloud, and the initial velocities of the matter shells, can evolve into either of 
these outcomes depending on the evolution chosen. In other words, there are permissible 
classes of evolutions, subject to the energy conditions and other conditions ensuring the 
physical reasonability, which give rise to either a black hole as the final product of collapse, 
or a naked singularity. 

It appears that the scenario presented above is repeated in nearly all known ex- 
act physically reasonable solutions of the field equations describing spherically symmetric 
gravitational collapse. In the known classes, such as radiation collapse [5], dust [11], per- 
fect fluids and more general forms of matter [7], and massless scalar field collapse [2], 
singularities both naked or covered do occur. Therefore, from the point of view of cosmic 
censorship, the important question currently is not the occurrence of naked singularities 
but whether these occurrences are permitted by the general relativity on a generic enough 
basis, and whether they are stable according to a suitably defined stability criteria. The 
issue of genericity is significant from the point of view that these phenomena could be 
impossibly rare in the space of solutions of the general theory of relativity. The known lit- 
erature does not indicate this clearly, because for all known important classes of solutions 
in the spherically symmetric gravitational collapse both the naked singularities and black 
holes occur for wide ranging situations. 

An important indicator in this direction, as pointed out earlier, is the imploding 
Vaidya model, where the singularity is naked for A < 1/8 and a black hole for A > 1/8, 
where the parameter A represents the initial data in the form of the rate of mass loss. 
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Thus, the singularities, both naked and covered, are stable against the pertubation of the 
parameter A, and the point A = 1/8 is the critical point indicating the transition from one 
phase to the other. A similar situation also occurs in the case of dust collapse where the 
pertubation in the initial density or velocity distribution within a certain domain does not 
alter the nature of the singularity. The analysis here in general has a significance in that 
the nature of the singularity is seen to be stable in a certain sense against the pertubation 
of the initial data. 

In this sense, we have characterized here wide new families of black hole solutions 
forming in spherically symmetric gravitational collapse, without using the assumption of 
the cosmic censorship hypothesis, but in terms of the possible evolutions of the initial data 
for the collapsing object. What we still do not know is the actual measure of each of these 
classes in the space of all possible evolutions which are allowed from a given general and 
arbitrary but physically reasonable initial data set. This is a problem related closely to the 
issue of stability of naked singularities. As is well-known, the stability in general relativity 
is a complicated issue because there is no well-defined formulation or specific criteria to test 
for stability. Fast evolving numerical codes for core collapse models may possibly provide 
further insights into this aspect. All the same, these classes appear to be generically arising 
in the collapse models considered here, at least within spherical symmetry, in that they 
are not an isolated phenomena but belong to a general family. Because, given any density 
and pressure profiles for the cloud, there exists an evolution which will lead to either a 
black hole or a naked singularity as desired as the end product of collapse. In this sense, 
both black holes and naked singularities do seem to arise generically as the end product for 
spherically symmetric gravitational collapse. Given the complexity of the field equations, 
if a phenomena occurs so widely in spherical symmetry, it is not unlikely that the same 
would be repeated in more general situations. 

It thus appears from the above considerations that the occurrence of singularities, 
naked or otherwise, is inherent in the theory of general relativity, and a distinction be- 
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tween these cases may not be possible through general relativity alone. What is essential 
is to examine how the perturbations and departures from sphericity would alter these con- 
clusions. Some efforts have been made to examine this issue [12]. It is also possible that 
the quantum effects near the naked singularities may help to preserve some kind of a quan- 
tum cosmic censorship, or these quntum effects could give rise to interesting signatures for 
naked singularities [13]. 
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